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Abstract 

We study the nonlinear-damping continuation of singular solutions of the critical and supercrit- 
ical NLS. Our simulations suggest that for generic initial conditions that lead to collapse in the 
undamped NLS, the solution of the weakly-damped NLS 

iV't(t,x) + AV' + iv^r" V + i'5iV'r"V = o, o < ^ < i, 



Ph \ is highly asymmetric with respect to the singularity time, and the post-collapse defocusing velocity 

■^ ■ of the singular core goes to infinity as the damping coefficient 5 goes to zero. In the special case 

,J^ . of the minimal-power blowup solutions of the critical NLS, the continuation is a minimal-power 

C^ \ solution with a higher (but finite) defocusing velocity, whose magnitude increases monotonically 

with the nonlinear damping exponent q. 

1. Introduction 

OO . The nonlinear Schrodinger equation (NLS) 



iVi(i,x) + AV + |Vr^V = 0, Vo(0,x)=^o(x)Gi/\ (1) 



where x = (xi, ..., x^) € M and A = dx^xi + • • 'dx^xa^ is one of the canonical nonlinear equations in 
physics, arising in various fields such as nonlinear optics, plasma physics, Bose-Einstein condensates 
(BEC), and surface waves. When (p — l)d < 4, the NLS is called subcritical. In that case, 
all H^ solutions exist globally. In contrast, both the critical NLS {p— l)d = 4 and the supercritical 
NLS {p — l)d > 4 admit singular solutions. Since physical quantities do not become singular, this 
C^ I implies that some of the terms that were neglected in the derivation of the NLS, become important 

near the singularity. 

The continuation of NLS solutions beyond the singularity has been an open question for many 
years. In 1992, Merle [1] presented a continuation of the explicit blowup solutions ^expiicit,o of 
the critical NLS, see (9), which is based on slightly reducing the power (L^ norm) of the initial 
condition. This continuation has two key properties: 

1. Property 1: The solution is symmetric with respect to the singularity time T^. 

2. Property 2: After the singularity, the solution can only be determined up to multiplication 

by a constant phase term e . 
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More recently, Merle, Raphael and Szeftel [2] generalized this continuation result to Bourgain- 
Wang singular solutions [3]. Note, however, that both the explicit solutions V'oxpiicit,o and the 
Bourgain-Wang solutions are unstable. 

In [4], Merle presented a different continuation, which is based on the addition of nonlinear 
saturation. Merle showed that, generically, as the nonlinear saturation coefficient goes to zero, the 
limiting solution beyond Tc can be decomposed into two components: A (^-function singular core 
that extends for Tc <t < T^, and a regular component elsewhere. 

In [5], Tao proved the global existence and uniqueness in the semi Strichartz class for solutions 
of the critical NLS. Intuitively, these solutions are formed by solving the equation in the Strichartz 
class whenever possible, and deleting any power that escapes to spatial or frequency infinity when 
the solution leaves the Strichartz class. These solutions, however, do not depend continuously on 
the initial conditions, and are thus not a well-posed class of solutions. Recently, Stinis [6] studied 
numerically the continuation of singular NLS solutions using the t-model approach. 

In [7] we analyzed asymptotically and numerically four potential continuations of singular NLS 
solutions: 1) a sub-threshold power continuation, 2) a shrinking-hole continuation for ring-type so- 
lutions, 3) a vanishing nonlinear-damping continuation, and 4) a complex Ginzburg-Landau (CGL) 
continuation. Our main findings were as follows: 

1. The non- uniqueness of the phase of the singular core beyond the singularity (Property 2) is 
a universal feature of NLS continuations. 

2. The symmetry with respect to the singularity time (Property 1) holds if the continuation 
model is time reversible and if it leads to a point singularity (i.e., if it defocuses for t > Tc). 
Therefore, it is a non-generic feature. 

Recently, the post-collapse loss-of-phase phenomena was demonstrated experimentally for intense 
laser beams propagating in water [8]. 

In this paper we further study the effect of small nonlinear-damping in the NLS 

iV'i(t,x) + AV' + |'0|P-V + ^5|V'|''"Vj = O, 0<5^l. (2) 

The addition of small nonlinear-damping is physical. Indeed, in nonlinear optics, experiments 
suggest that arrest of collapse is related to plasma formation, and nonlinear damping is used 
as phenomenological model for multi-photon absorption by plasma. In BEC, a quintic nonlinear 
damping term corresponds to losses from condensate due to three-body inelastic recombinations [9]. 
In addition, the nonlinear-damping term appears in the complex-Ginzburg-Landau (CGL) equation, 
which arises in a models of chemical turbulence, Poiseuille flow, Rayleigh-Bernard convection, 
Taylor-Couette flow, and superconductivity. 

In [7] we analyzed the continuation of the critical NLS with a vanishing critical nonlinear 
damping, i.e., equation (2) with p = q = 1 + 4/d. Since the NLS (2) is not time reversible, its 
solutions are asymmetric with respect to the time T^J^st ^^ which the collapse is arrested. In 
particular, in the limit 6n — )• 0+, the continuation of ipcxphcit,a{'t,r) is e*^Vexpiicit reo(^-^c — t,r), 
where n ~ 1.614. Hence, the defocusing velocity ko is higher then the focusing velocity a. When 
the initial condition leads to a loglog collapse in the undamped critical NLS, asymptotic analysis and 
numerical simulations suggest that the singular core expands beyond the singularity at a velocity 
that goes to infinity as 5 ^ 0+. 

The question that we address in this study is whether and how the results of [7] for q = p = 
1 + 4:/d will change in the following cases: 

1. The critical NLS with a supercritical damping exponent (i.e., q > p = 1 + 4:/d). 



2. The supercritical NLS with q > p > 1 + 4/(i. 

The paper is organized as follows. In Section 2 we provide a short review of NLS theory. In 
Section 3 we review previous rigorous, asymptotic, and numerical results on the effect of damping 
in the NLS. In Section 4 we show numerically that in the supercritical NLS, the nonlinear damping 
exponent q has to be strictly higher than the nonlinearity exponent p, in order to arrest the 
collapse. This is different from the critical case, where collapse is arrested for q> p- In Section 5 
we show that solutions of the supercritical NLS with a small nonlinear damping are asymmetric 
with respect to the arrest-of-collapse time Tj^rj-gsf ^^*^ that the post-collapse defocusing velocity 
of the singular core goes to infinity as the damping coefficient 8 goes to zero. In Section 6 we 
obtain similar results for the critical NLS with generic initial conditions that lead to a loglog 
collapse. In the special case of the minimal-power explicit blowup solution '<pc^p\idt,a{t-,'r) of the 
critical NLS, however, the continuation beyond the singularity is also defined for q < p, and is given 
by e'^Vex licit k( W^^'^'^ ~ ^j'^); where K{q) increases monotonically with q. Final remarks are given 
in Section 7. 

Overall, the qualitative effect of small nonlinear damping on the collapse is the same in the 
critical and supercritical NLS. One difference is that in the critical case collapse is arrested for q>p, 
whereas in the supercritical case collapse is only arrested for q > p. Another difference is that the 
distance between the damped solution around Tj^i-rest ^^^ the asymptotic profile of the undamped 
NLS is small in the critical case, but large in the supercritical case. Surprisingly, in the latter case, 
the profile near Tarrest appears to be given by a rescaled supercritical standing wave. 

2. Reviev^ of NLS theory 

The NLS (1) has two important conservation laws: Power conservation^ 



P(t) = P(0), P{t) = jl^l^dx, 



and Hamiltonian conservation 



of 



Hit) = HiO), H{t) = I iVV^I^dx - ^ / IV-r+^rfx. (3) 

The NLS (1) admits the waveguide solutions ijj = e^^R{r), where r = |x|, and R is the solution 

R"(r) + ^^—^R' -R + RP = 0, R'{0) = 0, P(oo) = 0. (4) 

r 

When d = 1, the solution of (4) is unique, and is given by 

i^W=(4^)"'""cosh-/<-')(^.). (5) 

When d > 2, equation (4) admits an infinite number of solutions. The solution with the minimal 
power, which we denote by R^^' , is unique, and is called the ground state. 



^We call the L^ norm the power, since in optics it corresponds to the beam's power. 



2.1. Critical NLS 

In the critical case [p — l)d = 4, equation (1) can be rewritten as 

iVt(t,x) + AV^ + |V|^/''V' = 0, Vo(0,x) = Vo(x)Gi/\ (6) 

and equation (4) can be rewritten as 

E!'{r) + '^—^E!-R + R^I^+^ = Q, R'{0) = 0, R{oo) = 0. (7) 

r 

Theorem 1 (Weinstein [10]). A sufficient condition for global existence in the critical NLS (6) 
is \\ipo\\2 ^ -fcr; whcrc Per = ||i?^'^^|||, and R^^' is the ground state of equation (7). 

The critical NLS (6) admits the explicit solution 

1 „/m / r \ ,^^,Lj.rl 

L 4 



^..„.,.(t.) = ^fl'»'(^)e'-- 



where 

m = n-t. .(„ = ^'-^., = ^^ (8b) 

More generally, applying the dilation transformation with X = a and the temporal translation 
Tc — > a^Tc shows that the critical NLS (6) admits the explicit solutions 

^„,„,„(,.) = _^i,(0) (_i_) .--.^-f , (9.) 

where 

L^{t) = a{T,-t), T^{t)= f -^ds = \-^ a>0. (9b) 

The explicit solutions (8)-(9) become singular at t = Tc. These solutions are unstable, however, 
as the have exactly the critical power for collapse. Therefore, any infinitesimal perturbation which 
decreases their power, will arrest the collapse. 

When a solution of the critical NLS, whose power is slightly above Pa, undergoes a stable 
collapse, it splits into two components: A collapsing core that approaches the universal ipji profile 
and blows up at the loglog law rate, and a non-collapsing tail {(f)) that does not participate in the 
collapse process: 

Theorem 2 (Merle and Raphael [11], [12], [13], [14], [15], [16], [17]). Let d = 1,2,3,4,5, and let ip 
he a solution of the critical NLS (6) that becomes singular at Tc- Then, there exists a universal 
constant m* > 0, which depends only on the dimension, such that for any ipQ S H^ such that 

V \m\\2 J 

the following hold: 

1. There exist parameters (T(t),xo(i),i(t)) G M x M'^ x M+, and a function 7^ G L^, such 
that 

'(/'(t,x) - Vij(t,x-xo(t)) — ^<^(x), t — > Tc, 



where 



«.,.) = ^R<»)(Ji) .-<•). (10) 



L'^/2(t) 

and R^^' is the ground state of equation (7). 
2. Ast — > Tc, 

^<" ~ ^ ( logiiS(r'-t)i ) ('"^'"^ "''^>- <"> 

^.^. Supercritical NLS 

In contrast to the extensive theory on singularity formation in the critical NLS, much less is 
known about the supercritical case 

#t(t,x) + AV' + |'(/'|P-V = 0, (p-l)d>4. (12) 

Numerical simulations and formal calculations (see [18, Chapter 7] and the references therein), 

and recent rigorous analysis in the slightly-supercritical regime < ^ 2 <C 1 [19] show that 

peak- type singular solutions of the supercritical NLS (12) collapse with a self-similar asymptotic 
profile ■0(3) where 

1 „ , . ir+ii^r-^ r /■* ds 



The blowup rate of L{t) is a square root, i.e., 



L{t) ~ K^T, - i, t -^ Tc, (14) 

where k > 0. In addition, the self-similar profile Q is the zero-Hamiltonian, monotonically- 
decreasing solution of 

3. Effect of linear and nonlinear damping - review 

In [20] , Fibich studied asymptotically and numerically the effect of damping on blowup in the 
critical NLS. He showed that when the damping is linear, i.e., 

iVt(i,x) + A'(/' + |V'|^/'^V' + «<5V' = 0, V(0,x) =Vo(x), (16) 

if the initial condition V'o(x) is such that the solution of (16) becomes singular for 6 = 0, then the 
solution of (16) exists globally only if 6 is above a threshold value 6c > (which depends on 0o)- 
Therefore, linear damping cannot play the role of "viscosity" in continuations of solutions of the 
NLS. When, however, the damping exponent is critical or supercritical, i.e.. 



z 



V'i(i,x) + AVj + (l + i<5)|V'|4/'^V = 0, 0<(5<C1, (17) 



or 



i^t(t, x) + AV' + \^p\^^'^^p + i^jfAl^- V = 0, < 6 <^ 1 q-l> 4/d, (18) 



respectively, then regardless of how small 5 is, collapse is always arrested. Therefore, Fibich sug- 
gested that nonlinear damping can "play the role of viscosity" in defining weak NLS solutions, i.e., 
we can define the continuation 

ih := lim 'ib^^\ (19) 

where ip^' is the solution of (17) or (18). 

Passot, Sulem and Sulem proved that high-order nonlinear damping always prevents collapse 
for d = 2. Antonelli and Sparber extended this result to d = 1 and d = 3: 

Theorem 3 ([21, 22]). The d- dimensional cubic NLS with nonlinear damping 

iV't + AV' + A|V'|V + «<5|V'l''"V = 0, AgM, 6>0, (20) 

where V'o(x) G H^{M.'^), 3 < q < oo if d = 1,2, and 3 < q < 5 if d = 3, has a unique global in-time 
solution. 

This rigorously shows that high-order nonlinear damping can play the role of " viscosity" . More 
recently, Antonelli and Sparber proved global existence for the case where the damping exponent 
is equal to that of the nonlinearity: 

Theorem 4 ([22]). Consider the cubic nonlinear NLS with a cubic nonlinear damping 

iV'i(t,x) + AV'+(l + i<5)|V'|V = 0, (21) 

where V'o(x) G H^{W^), x^o ^ L'^{W^), and d < 3. Then, for any 6 > 1, equation (21) has a unique 
global in-time solution. 

Theorem 4 does not show that critical nonlinear damping can play the role of viscosity. We 
note, however, that the asymptotic analysis and simulations of [7, 20] strongly suggest that the 
solution of (17) exists globally for any < 5 ^ 1. 

3.1. Explicit continuation oZ-^expiicit 

In [7], Fibich and Klein calculated explicitly the vanishing nonlinear-damping limit (19) of the 
explicit solution ■(/'explicit: 

Continuation Result 1 ([7]). Let ip^ '{t,r) be the solution of the NLS (17) with the initial con- 
dition 

i'oir) = ■(/'explicit (0, r ) , (22) 

see (8) . Then, for any ^ G M, there exists a sequence 6n — )• 0-|- (depending on 6), such that 

.i^+ ^ ^*' "^ - 1 V^:xpiicit,.(2Tc - 1, r)e^^ n<t<oo, ^^3) 

where ■i/'cxpiicit,fc is given by (9a) with a = k, 

K = 7r [Bi{0)A'i{s*) - AiiO)Bi{s*)] « 1.614, (24) 

Ai{s) and Bi{s) are the Airy and B airy functions, respectively, and s* ~ —2.6663 is the first negative 
root ofG{s) = V3Ai{s) - Bi{s). 



In particular, the limiting width of the solution is given by 

lim LW = |VV, l-'^^'- (25) 

,5^0+ ^ ' { K{t-Tc) Tc <t <cc, ^ ^ 

Therefore, the continuation is also an exphcit minimal-power solution, but with a higher (defo- 
cusing) velocity. In addition, the solution beyond Tc is only determined up to multiplication by an 
unknown constant phase e . 

3.2. Continuation for loglog collapse 

In [7] , Fibich and Klein showed asymptotically and numerically that the continuation of solutions 
that undergo a loglog collapse has an infinite-velocity expanding core, that is determined up to a 
multiplicative constant phase e : 

Continuation Result 2 ([7]). Let ^^{r) be a radial initial condition, such that the corresponding 
solution tp of the undamped critical NLS (6) collapses with the V'fi profile at the loglog law blowup 
rate at Tc- Let ip^ ' be the solution of the damped NLS (17) with the same initial condition. Then, 

lim ib^^'^ =ih, 0<t<Tc. 

5-5>0+ 

In addition, for any < 5 <^ 1, there exists 6 (5) £ M, and a function (p £ L'^ , such that 



lim 

<5-5>0+ 



'V(') (t, r) - V|j(2r, - t, r; 5)e''^'^] ^ ^{r), t -^ T,+, 



where ipR is given by (10) with some function L{t;5), such that 

lim lim L(t;5) = 0, lim lim Lt(t;5) = oo, lim 9(6) = oo. 

Remark 1. We use the terminology Continuation Result, in order to emphasize that the proofs of 
Continuation Results 1 and 2 are based on asymptotic analysis and numerical simulations, and are 
not rigorous. 

4. The critical damping exponent 

The rigorous theorems in Section 3 suggest that solutions of the NLS (2) with ^o(x) G H^ 
always exist globally when p < q. These theorems, however, do not cover the case p = q. The 
asymptotic analysis and simulations of [7, 20] strongly suggest that in the critical case, when p = q 
the solution always exists globally, see Continuation Results 1 and 2. Since there are no rigorous 
and asymptotic results for the supercritical case with p = q, we study this case numerically. 

Consider the one-dimensional NLS with p = q 

iMt, x) + ^^^ + (1 + i5)| Vr" V = 0, (26a) 

with the perturbed solitary-wave initial condition 

V'o(x) = 1.05Rp{x), (26b) 



where Rp{x) is given by (5). In Figure 1 we solve (26) with 5 = 5 • 10 ^, and plot 



m 



^(0,0) 



V^(t,o) 



(p-l)/2 



(27) 



In the critical case p = 5, the collapse is arrested after focusing by ~ 10. In the supercritical 
case p = 7, however, the collapse is not arrested after focusing by 10^. This and similar simulations 
suggest that unlike the critical case, in the supercritical case, the condition p < q is necessary for 
ensuring global existence in (2). 




Figure 1: Solution of (26) for 5 = 5 • 10 with p — 5 (dashes) and p = 7 (sohd). 



5. Supercritical NLS 

We now consider the effect of small nonlinear damping in the supercritical NLS. Let t(j^ '{t,x) 
be the solution of the one-dimensional supercritical damped NLS {d = 1, p = 7 , q = 9) 



iMt, x) + Vxx + 1-01 V + i<5|V'l V = 0, 



with the initial condition 
Let 



ipo{x) = 1.3e 



^(5) 



max - arg max | 

t 



',{S) 



it,x) 



(28a) 

(28b) 
(29) 



denote the time at which the focusing is maximal. 

In Figure 2 we plot the solution of (28) for various values of 6. In all cases, the collapse is arrested 
in a highly asymmetric way with respect to T!^„csf ^^ addition, the post-collapse defocusing rate 
of the singular core "appears" to increase to infinity as (5 — )• 0-|-. 

In Figure 3 we compare the profile of the solution of (28) with 6 = 10"^ with the supercritical ^q 
and ipR profiles, where 



\'<pR{t,x)\ 



1 



-R 



f X 



Lf^-'\t) \Ln{t)J' 



LR{t) 



RiO) 



V'(^)(t,o) 



(p-l)/2 



(30a) 




Figure 2: Solution of (28) for 5 = (solid), 5 = 5- 10"^ (dashes-dots), 5 = 7.5 • 10"^ (dots), and 5 = 10"^ (dashes). 



\'ipQ{t,x)\ 



Lf'-'\t) 



Q 



X 



Lgit) 



Lgit) 



Q{0) 



Tp(^){t,0) 



(p-l)/2 



(30b) 



and R and Q are the solutions of (4) and (15), respectively, with d = 1 and p = 7. The NLS solution 
initially approaches the ijjq profile, see Figure 3(a-c). This is to be expected, since when 6 = the 
solution collapses with the ^q profile, see Section 2.2. As the solution approaches T^Jcsf however, 
the collapsing core moves away from ipQ and towards ipji, see Figure 3(d), and it remains close 
to ipji for a "short time" after T^i-rest' ^^^ Figure 3(e). Eventually, as the collapsing core continues 
to defocus, it interacts with its tail and "loses" its ^r profile, see Figure 3(f). 

Next, we repeat the above simulation with a higher nonlinear damping exponent {q = 11). 
Specifically, we solve the NLS 



i^Ptit, x) + V'xx + IV'IV + i'JiV'l > = 



(31a) 



with the initial condition 



ipoix) = l-3e 



(31b) 



Figure 4 and 5 show that the qualitative behavior of the solution is exactly the same as that of the 
solution of (28). 

Therefore, we conclude that solutions of the supercritical NLS (2) with q > p > 1 + 4/d 
and < (5 < 1: 



.(<5) 



1. Exist globally. 

2. Are highly asymmetric with respect to T^^^rrcst- 

3. The post-collapse velocity of the defocusing core goes to infinity as 6 

4. The asymptotic profile around T!^„cst is ipR, and not ^q. 

.(-5) 



0+. 



The fact that as i — t- T.^j.^^s^^ the profile is not given by ipg is not surprising, since the nonlinear 

damping perturbation obviously has a significant effect near T^^j.^^^.^^, and therefore there is no 
reason why it should not change the solution profile. What is surprising is that the profile 
changes to the supercritical ipR profile. As far as we know, this is the first observation in 
which the asymptotic profile in the supercritical NLS is given by the supercritical ipR profile.^ 



^The standing-ring solutions of the undamped supercritical NLS with p = 5 and d > 1 also collapse with the 
tpR profile [23, 24, 25, 26]. In that case, however, tpR is the asymptotic profile of the critical one-dimensional quintic 

NLS. 




lOL 



lOL 



0.25 



0.2375 



0.2377 



Figure 3: Solution of (28) for 5 = 5 • 10^^ (solid), and the fitted |V'q| (dots) and IV-fll (dashes), (a) t « 0.2, L « 
0.38 (b) i ^ 0.227, L ^ 0.2 (c) t ^ 0.2368, L ^ 0.05 (d) t ^ 0.2376, L ^ 0.0118 (e) t « 0.23764, L « 0.0119 
(f) t ~ 0.23767, L ~ 0.135. The stars in (g) and (h) denote the values of t and L{t) for the data in subplots (a)-(f). 




Figure 4: Solution of (31) for 5 = (solid), S = 5- 10"* (dots), and S = 10"^ (dashes). 



6. Critical NLS 

6.1. Continuation of loglog collapse 

In [7], we studied the effect of nonlinear damping in the critical NLS with p = q with initial 
conditions that lead to a loglog collapse, see Continuation Result 2. We now consider the case p > q. 
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0.2226 



Figure 5: Same as Figure 3 for the solution of (31). (a) t ^ 0.19, L ^ 0.33 (b) i ^ 0.211, L ^ 0.22 (c) t ^ 0.221, L : 
0.033 (d) t « 0.222, L ^ 0.02341 (e) t f» 0.2222, L ^ 0.02343 (f) t f» 0.2223, L ^ 0.0267. 



Consider the damped one-dimensional critical NLS {d = 1, p = 5, q = 7) 

iMt, x) + Vxx + 1-01 V + i6\^f^ = 0, 



with the initial condition 



ipo{x) = l-6e 



(32a) 



(32b) 



whose power is 4% above the critical power for collapse. When 6 = 0, the NLS solution collapses 
with the ipR profile at the loglog blowup rate. 

In Figure 6 we solve (32) for various values of 5. In all cases, the collapse is arrested in a highly 
asymmetric way with respect to T!^„esf ^^^ addition, the post-collapse defocusing rate appears to 
increase to infinity as 5 — t- 0+. This qualitative behavior is as in the case p = q, see Continuation 
Result 2. Therefore, we conclude that the qualitative behavior for q = p and for q > p is the same. 

In Figure 7 we compare the profile of the solution of (32) with 6 = 10~^ with the best-fitting 
critical ipR profile, see (30a). The NLS solution initially approaches the ipji profile, see Figure 7(a- 
c). This is to be expected, since when 6 = the solution collapses with the ipR profile, see 
Theorem 2. As the solution approaches T!^„cst^ however, the collapsing core moves away from tpji, 
see Figure 7(d-e). Unlike the supercritical case, however, the solution profile near T!^„cst ^^ still 
"close" to ipR. This is because in the critical case, perturbations arrest the collapse when they are 
still small compared with the nonlinearity and diffraction [27]. Eventually, as the collapsing core 
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L(t) 




0.2 0.22 0.24 0.26 0.28 0.3 
t 



Figure 6: Solution of (32) for 5 = 10"'' (solid), S = 2.5 ■ 10"" (dots) , and 5 = 5 ■ 10"* (dashes). 



continues to defocus, it interacts with its tail and "loses" its tpR profile, see Figure 7(f). 
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Figure 7: Solution of (32) for S = 10"^ (solid), and the fitted I^jr] (dashes), (a) t ?a 0.205, L ^ 0.45 (b) t ?a 0.26, L : 
0.24 (c) t « 0.28514, L ^ 0.0037 (d) t ^ 0.285146, L ?» 0.0032 (e) t ^ 0.285147, L ^ 0.0032 (f) t ^ 0.28515, L ; 
0.0046. The stars in (g) and (h) denote the values of i and L{t) for the data in subplots (a)-(f). 



In summary, nonlinearly-damped loglog solutions of the critical NLS with q > p have the 
following properties: 
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1. The solutions are highly asymmetric with respect to ^arrest- 

2. The post-collapse defocusing velocity goes to infinity as (5 — )• 0+. 

3. The asymptotic profile near T^^^^^t i^ "slightly" different form ipR 

6. 2. Continuation of '(/'explicit 

Consider the critical NLS with nonlinear damping 



i^t, x) + AV' + IV'I^/'^^ + iJI^I^- V = 0, 
and the initial condition 

V'o(?') = '0explicit(O,r). 



0<(5«1, 



(33a) 



(33b) 



When 5 = 0, the solution is given by ^explicit, see equation (8). In [7], we calculated explicitly the 
continuation of V'expUcit when q = p, see Continuation Result 1. We now consider the continuation 
for q 7^ p.^ 

As in [7], we can use modulation theory [27] to approximate equation (33) with a reduced system 
of ordinary-differential equations. 



Lemma 1. Let 



m 



^^(0,0) 



V^(t,o) 

where ip is the solution of equation (33). Then, as 5 - 
the reduced equations 



2/d 



2Cq5 



1 



0-|-, the evolution of L{t) is governed by 

m 



M L(5-iW2' 



Ltt{t) 



L3 



subject to the initial conditions 

m = 0, 

where 

v{l3) = 



L(0) = 0, Li(0) 



(34a) 
(34b) 



0, 



/3>0, 
/3<0, 



2A| 

M ' 



M = - 



^li^W 



Ar= lim eM'^-^)/2^(°)(r 



V-idr, c, = \\R^'XXl 



and R^^' is the ground state of (7). 

Proof. In [20] it was shown that the reduced equations for the damped NLS (33a) are given by 

u{l3) 2cg6 1 ^ _ m 



fit{t) 



L2 



M L(l'l)d/2' 



Lt 



L3 ■ 



(35) 



In addition, the initial conditions for the reduced equations (35) that correspond to the initial 
condition (33b) are /3(0) = 0, L(0) = 0, and Lt(0) = 0, see [7]. Since /3(0) = 0, and since ft < 0, 
then /3(t) < 0. Hence, ;/(/?) = 0. Therefore, the reduced equations are given by (34). D 



^Since -^explicit has exactly the critical power for collapse, any amount of damping will arrest the collapse. There- 
fore, the continuation of '(/'explicit can also be defined for q < p. 
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The reduced equations variable L{t) is the sohition width, and is also inversely proportional 
to the solution amplitude. The reduced equations variable /3(i) is a measure of the acceleration 
of L{t), and is also linearly proportional to the excess power above Per of the collapsing core. 

Since modulation theory is not rigorous, in Figure 8 we compare the numerical solutions of 
the reduced equations (34) and the NLS (32). This comparison shows that the two solutions are 
in excellent agreement, thus providing a strong support to the validity of the reduced equations. 
Therefore, in what follows we study asymptotically and numerically the limit 6 — > 0+ within the 
framework of the reduced equations, which is considerably easier than studying the limit S — > 0+ 
of the nonlinear ly-damped NLS. 




Figure 8: Solution of the reduced equations (34) [solid], and of the NLS (33) [dashes], for S = 2.5 ■ 10 ^ and d = 1. 
The two curves are indistinguishable, (a) g = 1. (b) q — 3. (c) q — 5. (d) q = 7. 



The extension of Continuation Result 1 to (7 7^ p is as follows. 

Continuation Result 3. Let ^l)^^\t,r) be the solution of the NLS (33). Then, for any G 
there exists a sequence 6n — > 0+ (depending on 9), such that 



lim ^(''")(t,r) 



V'explicit(i,'r) 

V':xpiicit,K(q)(2rc-t,r)e^^ 



< t < Tc, 

Tr < t < 00. 



(36) 



In particular, the limiting width of the solution is given by 



lim L(t) 



Tc-t 
K{q){t-Tc 



0<t<Tc, 
Tc<t< 00. 



(37) 



Proof. We only provide an informal proof, using the reduced equations (34). As (5 — t- 0+, 

/3t(t) — > 0, see equation (34a). Therefore, since /3(0) = 0, then /3(t) — ;• 0. Hence, Lu{t) — )■ 0. 
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Therefore, lim5_j.o+ L{t) is linear in t. Since lim5_j>o+ L{Tc) = 0, it follows that 



lim L(t:6) 



t, 



t<Tc, 

n<t. 



K{q){t-T,), 
The loss of phase follows from the fact that liuit^Tc arg('0cxpiicit(i)O)) = oo ,see [7]. D 

The result of Continuation Result 3 can be explained as follows. By continuity, lini5_>>o+ tp^ ' = 
V'expiicit for < t < Tc- The limiting solution for t > Tc is an NLS solution that becomes singular 
as i — )• Tc+, and has exactly the critical power at the singularity. Hence, the limiting solution is a 
minimal-power solution. Therefore, it has to be given by V'expiicit,^ [li 4]. 

In Figure 9 we solve the reduced equations (34) with 6 = 10~^, and observe that: 

1. The limiting solutions are indeed linear for t < Tc and t > Tc- 

2. The continuation is asymmetric with respect to Tc- 

3. The post-collapse slope K{q) increases with q. In [7] we showed that the jump discontinuity 
in lim5^o+ -^t ^-t Tc is related to the increase of the Hamiltonian as the limiting solution 
passes through the singularity. As q increases, damping affects become more pronounced, 
hence there is a larger increase of the Hamiltonian, hence of the post-collapse slope. 

4. When q = 1, K{q =!) = !, i.e., L{t) is symmetric with respect to Tc- Therefore, the linear 
damping continuation of ^explicit is symmetric with respect to Tc, even though the problem is 
not time-reversible. 



Note that the value of n^q 
Result 1. 



1 + 4/d) « 1.614 was computed analytically in Continuation 




Figure 9: Solution of the reduced equations (34) with 5 = 10 ^ and various values of q. (a) L{t). (b) K.{q). 



7. Final remarks 

In this study we used numerical simulations to study the effect of small nonlinear-damping 
on singular NLS solutions. These simulations suggest that the effect of small nonlinear damping 
is qualitatively the same in the critical NLS with generic initial conditions that lead to a loglog 
collapse with the ipR profile, and in the supercritical NLS with generic initial conditions that lead 
to a square-root collapse with the V'Q profile. Moreover, the qualitative effect of nonlinear damping 
is independent of the value of q, so long as q > p in the supercritical case and q > p in the critical 
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case. Thus, because nonlinear damping destroys the NLS time reversibihty, the nonhnearly-damped 
solution is highly asymmetric with respect to the arrest-of-collapse time T^Jcst- '^^^ post-collapse 
defocusing velocity Lt{t) of the singular core goes to infinity as 6 — t- 0+, since the focusing velocity 
before the singularity goes to infinity for loglog and square-root blowup rates, and since nonlinear 
damping increases the Hamiltonian,^ hence the "kinetic energy". 

Around ^arrest ' *^^ Collapsing core of the singular core moves away from the asymptotic profile 
of the undamped solution. In the supercritical case the difference between the solution profile 
and ipQ for t w T^rLst is large. This is intuitive, since damping effects have a large effect when they 
arrest the collapse. In the critical case, however, the difference between the solution profile and ipn 
for t ~ ^arrest ^^ minor. This is because critical collapse has the unique property that it can be 
arrested by small perturbations [27]. 

Surprisingly, in the supercritical case the profile of the nonlinearly- damped solution near T!^„Qst 
appears to be given by the supercritical tp^ profile. To the best of our knowledge, this is the first 
observation of a solution of the supercritical NLS that approaches the supercritical ipR profile. 
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